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ABSTRACT 

We present new radio, optical, and X-ray observations of three Ultraluminous X-ray sources (ULXs) 
that are associated with large-scale nebulae. We report the discovery of a radio nebula associated 
with the ULX IC342 X-l using the Very Large Array (VLA). Complementary VLA observations of the 
nebula around Holmberg II X-l, and high-frequency Australia Telescope Compact Array (ATCA) and 
Very Large Telescope (VLT) spectroscopic observations of NGC5408 X-l are also presented. We study 
the morphology, ionization processes, and the energetics of the optical/radio nebulae of IC342 X-l, 
Holmberg II X-l and NGC5408 X-l. The energetics of the optical nebula of IC342 X-l is discussed in 
the framework of standard bubble theory. The total energy content of the optical nebula is 6 x 10 52 
erg. The minimum energy needed to supply the associated radio nebula is 9.2 x 10 50 erg. In addition, 
we detected an unresolved radio source at the location of IC342 X-l at VLA scales. However, our Very 
Long Baseline Interferometry (VLBI) observations using the European VLBI Network likely rule out 
the presence of any compact radio source at milli-arcsecond (mas) scales. Using a simultaneous Swift 
X-ray Telescope measurement, we estimate an upper limit on the mass of the black hole in IC342 
X-l using the "fundamental plane" of accreting black holes and obtain Mbh < (1-0 ± 0.3) x 10 3 M . 
Arguing that the nebula of IC342 X-l is possibly inflated by a jet, we estimate accretion rates and 
efficiencies for the jet of IC342 X-l and compare with sources like S26, SS433, IC10 X-l. 
Subject headings: black hole physics — accretion, accretion disks — ISM: bubbles, jets and outflows 
— X-rays: binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources are variable non-nuclear 
X-ray sources in external galaxies with luminosities 
greatly exceeding the Eddington luminosity of a stellar- 
mass compact object, assuming isotropic emission 
([Colbert fcMushotzky|[l999tlKaaret e"t~aT][200l . The ir- 
regular variability, observed on time scales from seconds 
to years, suggests that ULXs are binary systems contain- 
ing a compact obj ect that is either a stellar- mass black 
hole with beamed ([King et al.ll2001l:lKording et alll2Tj02ft 
or super-Eddington emission (Bcgehnasl |2002); or an 
intermediate- mass black hole (IMBH). IMBHs have been 
invoked in contexts ranging from the remnants of Pop- 
ulation III ([Madau fc R ces 2001) stars to the formation 
of su permassive black holes (SMBHs) (jPtak fc Griffiths 
1999); SMBHs may form thro ugh the hierarchical me rger 
of lower mass black holes (eg. lEbisuzaki. et al.ll2001ft . 

Several ULXs show emission-line optical nebulae, 
which can be used as a calorimeter to infer the to- 
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tal in t rinsic power of the U LX ([Pakull fc Mirionil [20021 
120031: iRussell et all 12011ft . In general, the nebu- 
lae around ULXs are either photoionized due to the 
high X-ray and UV lum i nosity of the compact object 
(IPakull fc Mirionil [2001 [20031 IKaaret. Ward fc ZezasI 
120041: iKaaret fc Corbel! [20091 iMoon et all 12011ft . or 
shock- io nized driven by j ets, o utflows, and/o r disk 
winds (IPakull fc Mirionl l2002t iRoberts et all 120031: 
lAbolmasov et a l. 2007). In several cases, two-component 
optical line profiles are present, indicating a mix- 
ture of the two mechanisms or alternatively the nar- 
row line could be due to the shock precursor. An- 
other common feature is the presence of the high- 
ionization He II emission line. It can have vari- 
ous origins: the nebula, the donor, the accretion 
disk, or a disk wind. The ionizing photon rate 
needed to produce the narrow high-excitation He n 
line of the nebulae in photoionized sources indicates 
that their X-ray emissi o n is a t most mildly beamed 
(IPakull fc Mirionil [2002[ [2001 IKaaret Ward fc ZezasI 
120041: IKaaret fc Corbelll2009l:lMoon et al.ll2011ft~ 

Only a handful of radio detections o f ULXs have 
been made so far including NGC5408 X-l jKaaret etall 
[20031: ISoria et al l I2006al ILang et al.1 [2007ft. Ho I I X-l 
dMiller et all 12005ft. and MF16 ([van Dvk et all [l99l 



ILacev. Duric fc Gossl [1 9971 , if not con sidered a super- 
nova remnant (Matonick & Fesenl 11997ft . Most of these 
sources show large nebulae (>~50 pc) that are likely 
powered by continuous energy input from the ULX. 
Shock-dominated ones are probably powered in the same 
manner as the W50 nebula is pow ered by the Galactic 
binary SS433 ([Dubner et all 11998ft . However, the ULX 
radio nebulae require greater total energy content than 
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W50. A similarly powerful n ebula, S26, w as found by 
iPakull Soria & Motchl (|2010[ ) in optical and iSoria et all 
( 20101) in radio. For other possib le radio associations 
with ULXs, we refer the r e ader to Sancb iez^Sutil et al. I 
(|200ll : ISoria et al.l (|2006bD : IMezcua fc Lobanovl (|2oToi r 

The X-ray spectra of ULXs share some similar 
properties with the canonical black hole states of 
Galactic black hole binaries (GBHBs) . A num- 
ber o f ULXs show state transi ti ons (iKubota et all 
2001: lFen7fc Kaaretl l2?M 12001 IGodet et al.l 1200ft 
Grise et al. I 120101 : iServillat et al.l 120111 ). InGBHBs 
during the X-ray hard state, the sources are associ- 
ated with self-absorbed compa ct jets ([Corbel et alll200~4l : 
iFender. Belloni fc~G allo 2004) . Given the similarities be- 
tween ULXs and GBHBs, it is interesting to investigate 
the presence of such compact jets for ULXs with hard 
X-ray spectra. 

In this paper, we present new radio and optical obser- 
vations of two ULXs that are associated with large-scale 
nebulae: Holmberg II X-l (Ho II X-l) and NGC5408. In 
addition, we present discovery of a radio nebula associ- 
ated with IC342 X-l. In Section 2, we describe obser- 
vations, data analysis and results. The energetics of the 
optical and radio nebulae and the jet properties of IC342 
X-l are investigated in Section 3. In Section 4, we briefly 
summarize our results. 

2. OBSERVATIONS, ANALYSIS, AND RESULTS 

2.1. VLA observations of IC342 X-l and Ho II X-l 

Observations of IC342 X-l and Ho II X-l were car- 
ried out using the C- and B-array configurations of the 
Very Large Array (VLA) of the National Radio Astron- 
omy Observatory (NRAO). The observations were made 
at 4.8 and 8.5 GHz (VLA program code: AL711) and 
the details are summarized in Table 1. Data calibration, 
combination in the (u,v) plane and imaging were car- 
ried out using the NRA O Astronomica l Image Process- 
ing System (AIPS: e.g., iGreisenl (|2003l) l We adjusted 
the Robust weighting parameter between and -2 to 
bring out the fine-scale radio emission and parameters 
are identified in image captions. 

2.1.1. Radio detection of IC342 X-l 

Figure Q] shows the 5 GHz VLA B- and C- array com- 
bined image of IC342 X-l wit h Robust=0 weighting over- 
laid on the Ha HST image of lFeng fc Kaaretl ()2008l i . Ex- 
tended radio emission is present surrounding the position 
of IC342 X-l. Diffuse emission is detected at up to the 
10-er level with peak intensity of ~120 /iJy beam -1 and 
the estimated total flux density in the nebula is ~ 2 mJy. 
The corresponding luminosity at a distance of 3.9 Mpc 
(jTikhonov fc GalazutdinaTOll20T0l ) is L neb = 1.8 x 10 35 
erg s _1 ; with L = vL v . Wc find that the size of the neb- 
ula is 16" x 8" , which corresponds to 300 pc x 150 pc. 
We estimate the volume of the nebula by taking a sphere 
with a diameter of 12" . Comparing the 5 GHz radio map 
to the Ha image, the optical and the radio nebulae are 
both elongated in the NE-SW direction, possibly exhibit- 
ing a shock- front. The size of the radi o nebula is consis- 
tent with the size of the optical nebula (|Pakull fc Mirionil 
I20i51 IGrise et al 1120061 iFeng h Kaaret]l2008D . 

Figure [2] illustrates a more uniformly weighted image 
of the 5 GHz radio emission surrounding IC342 X-l. We 



weighted the radio image with a robust parameter of -2 
in order to resolve out the diffuse nebular emission and 
show the fine scale structure. The strongest radio emis- 
sion appears towards the NE of the ULX and is coinci- 
dent with the strongest Ha emission. The radio emission 
extends farther to the NE in regions of little or no Ha 
emission. 

Figure [2] also reveals an unresolved radio source at the 
location of IC342 X-l. This unresolved source is de- 
tected with a flux density of ~ 96.3 /xJy at the 6.5-<r level. 
The corresponding luminosity is 8.8 x 10 33 erg s _1 . The 
HST position of the ULX is RA=03h45m55.61s, Decl.= 
+68°04'55. 3" (J2000.0) with the 90% positional errors of 
0.2 arcsec (|Feng fc Kaaretl 12008). We obtained the po- 
sition of the radio peak of RA=03h45m55.54s, Decl.= 
+68°04'55.18" using the maxfit task of AIPS. The dis- 
tance between the two positions is ~0.4 arcsec. Using the 
AIPS task jmfit, we estimate a radio positional uncer- 
tainty of ~0.1 arcsec at one sigma. Therefore the optical 
and radio positions are consistent using 90% positional 
errors. 

2.1.2. Multi-frequency radio observations of Ho II X-l 

The radio nebul a of Ho II X - l was first detected at 
1.4 and 5 GHz by IMiller etaLl (|2005l ). Here, we have 
conducted observations of Ho II X-l at 5 and 8 GHz 
in order to constrain the shape and spectrum of the 
radio nebula. Figure [3J shows the VLA images of 
Ho II X-l overlaid on HST H e II and H/3 images from 
iKaaret. Ward fc Zezai (|2004f ). On the left, the un- 
weighted (Robust=0) 5 GHz B-array image is shown. 
The right panel shows a Robust=— 1 image made at 
8 GHz using the C-array configuration. The asymmetric 
morphology of the nebula might indicate some outflows 
or ambient density gradient to the West. 

Previo usly, the radio sp ectrum was not well con- 
strained; Mi ller et all ([2005) had high uncertainties on 
the flux at 5 GHz of Ho II X-l. However, we now have 
flux density measurements at th ree frequenci e s (1.4 , 8 
and 5 GHz, using the image from IMiller et all (|2005h at 
1.4 GHz) and Figure |U illustrates the fitted three-point 
spectral index of a — —0.53 ± 0.07, (S ~ v a ). The radio 
spectrum is consistent with optically thin, synchrotron 
emission, further constraining the nebular nature of the 
radio counterpart of this ULX nebula. The integrated 
flux density of the nebula is 1.05±0.10 mJy, 613±61 /xJy 
and 395±40 ^Jy at 1.4, 5 and 8 GHz. 

2.2. Swift/XRT and EVN observations of IC342 X-l 
2.2.1. VLBI observation 

We conducted simultaneous X-ray and radio measure- 
ments of the VLA core (see Sec. 2.2.1.) to test its 
compactness in order to estimate the mass of the BH 
(see Sec. 3.4) with the minimum uncertainty. To check 
whether the VLA core is consistent with a compact jet 
at 10 mas scale, we carried out European VLBI Net- 
work (EVN) observations (EVN program code: EC032) 
at 1.6 GHz. The 12-hour observations were accommo- 
dated on June 15, 2011 from 03:30 to 15:30 (UT), simul- 
taneously with the Swift/XRT observations. The partici- 
pating VLBI stations were Effelsberg (Germany), Jodrell 
Bank Lovell Telescope, Cambridge (United Kingdom), 
Medicina, Noto (Italy), Torun (Poland), Onsala (Swe- 
den), Urumqi (P.R. China), Svetloe, Zelenchukskaya, 
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TABLE 1 

Summary of Observations 



Source 


Instrument 


Contig 


Central 


On-sourcc 


O bservation 


Flux 








Frequency 


Time 


Date 




1C342 X-l 


VLA 


B 


4.8 GHz 


3.2 h 


6 Dec 07 


2.0±0.1 mjy 


IC342 X-l 


VLA 


C 


4.8 GHz 


3.2 h 


25 Apr 08 


* 


Ho II X-l 


VLA 


B 


4.8 GHz 


3.2 h 


8 Dec 07 


613±61/iJy 


Ho II X-l 


VLA 


C 


8.5 GHz 


3.2 h 


21 Apr 08 


395±40^Jy 


NGO5408 X-l 


atoa 


6D 


5.5 GHz 


12 h 


23 Aug 09 


226±33AtJy 


NGC5408 X-l 


ATCA 


6D 


9 GHz 


12 h 


23 Aug 09 


137±36^Jy 


NGC5408 X-l 


ATCA 


6D 


17.9 GHz 


12 h 


23 Aug 09 


76±20/iJy 


NGO5408 X-l 


ESO VLT 




575-731 nm 


0.7 h 


12 Apr 10 




IC342 X-l 


EVN 




1.6 GHz 


12 h 


15 Jun 11 


< 21/xJy** 


IC342 X-l 


Swift/XRT 


PC 


0.3 - 10 keV 


9.4 ks 


15 Jun 11 


2.67*** 



Note. — *B and C configuration data were combined. **Three-sigma upper limit. ***Unabsorbed flux in units of 10 12 erg cm 2 s 1 . 




03 h 45 m 58 s .O 57 s . 56 s . 55 s . 54 s .O 53 3 .0 



Fig. 1. — VLA 5 GHz image of IC342 X-l overlaid on the H« HST image. Contours represent radio emission and are drawn at 3, 4, 5, 
6, 7, 8, 9, and 10 times the rms noise level of 11 fijy beam -1 . The peak brightness is 122.4 /xjy beam . The resolution of the image is 
2{'3 X 1'.'6 at PA=— 13° and the image was made with Robust=0 weighting. The sign '> <' marks the X-ray position of the ULX. 



TABLE 2 

Hubble images and source parameters 



Name 


Narrow Band Filter 


Centered on 


Ref. 


Distance 


Optical 


Radio 












diameter 


diameter 


10342 X-l 


F658N 


Ha 


1 


3.9 Mpc 


190 pc 


225 pc 


Ho II X-l 


FR462N 


He II 


2 


3.39 Mpc 


45 pc 


60 pc 


Ho II X-l 


FR505N 


H/3 


2 




101 pc 




NGC5408 X-l 


F502N 


[O m]A5007 


3 


4.8 Mpc 


60pc 


40 pc 



Note. — References: 1 -[Peng & Kaarct (2008), 2 - Kaarct. Ward & Zczas (20QD, 3 - IGrise et aT~l (I2012T) 
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Fig. 2. — VLA 5 GHz image of IC342 X-l overlaid on the Ha HST image shown in Figure [T] This image was made with Robust=-2 
weighting and therefore the largest extended features are missing. The resolution is 1'.'6 X l'.'l at PA=— 19°. Contours represent radio 
emission and are drawn at 3, 4, 5, and 6 times the rms noise level of 15 /xjy beam -1 . The peak brightness is 96.3 fijy beam -1 . The sign 
'> <' marks the X-ray position of the ULX. 




0.001 0.002 0.003 0.004 0.005 0.006 0.001 0.00S 0.01 0.015 0.02 0.02S 0.03 



Fig. 3. — Left: The 5 GHz VLA B-array image of Ho II X-l overlaid on the He II HST image. Contours represent radio emission at 
levels of 3, 4, 5, 6 and 7 times the rms noise level of 14 fijy beam -1 . The peak brightness is 114 fi.]y beam -1 . The image was made using 
Robust=0 weighting and the resolution is l'/5 X l'.'O at PA=34°. Right: The 8.5 GHz VLA C-array image of Ho II X-l overlaid on the 
H/3 HST image. Contours represent radio emission at levels of 3, 4, 5, 6 and 7 times the rms noise level of 15 /^Jy beam -1 . The peak 
brightness is 145 /xjy beam . The image was made using Robust=-l weighting and the resolution is 2'/36 X 1.75 at PA=3°. The sign 
'> <' marks the X-ray position of the ULX. The North direction is up in both images. 
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Fig. 4. — The radio spectrum of the nebula of Ho II X-l. The 

best fit spectral index is a = — 0.53±0.07 for a flux density, S oc v a . 

Badary (Russia) and the phased array of the Westerbork 
Synthesis Radio Telescope (WSRT; The Netherlands). 
The aggregate bitrate per station was 1024 Mbps. There 
were eight 8 MHz intermediate frequency channels (IFs) 
in both left and right circular polarisations. 

The source was observed in phase-reference mode. 
This allowed us to increase the coherent integration time 
spent on the target source and thus to improve the 
sensitivity of the observations. Phase- referencing in- 
volves regularly interleaving observations between the 
target source and a nearby, bright, and compact refer- 
ence source. The delay, delay rate, and phase solutions 
derived for the phase-reference calibrator (J0344+6827) 
was interpolated and applied to the target within the 
target-reference cycle time of 5 min. The target source 
was observed for 3.5-min intervals in each cycle. 

AIPS was used for the VLBI data c alibration following 
standard procedures (|Diamondl lT995) . The visibility am- 
plitudes were calibrated using system temperatures and 
antenna gains measured at the antennas. Fringe-fitting 
was performed for the calibrator sources using 3-min so- 
lution intervals. The calibrated visibility data were ex- 
porte d to the Caltech Difmap program (Shcperd ct al. 
119941 ) used to make a naturally weighted VLBI image. 
The achieved 1-a rms noise level was 7 /iJy beam -1 and 
we did not detect any source above the 3-cr noise level 
in a field of view of 700 x 700 mas. Therefore, the VLA 
core is likely to be a clump of emission from the nebula. 
Our EVN observation places a 3-<r upper limit on the 
flux density of F v < 21 fiJy and on the luminosity of the 
putative compact jet of < 1.9 x 10 33 erg s — 1 . 

2.2.2. Swift observation 

The Swift X-ray Telescope (XRT) obtained 9394 sec- 
onds of good exposure in its photon-counting (PC) mode 
on June 15, 2011, from 02:57:31 to 17:47:47 (UT). We 
retrieved level two event files and used the default data 
screening methods as described in the XRT user's guide0. 
We extracted an X-ray spectrum for the source using a 
circular extraction region with a radius of 20 pixels (cor- 



responding to 90% of the PSF at 1.5 keV); background 
was estimated from a nearby circular region with a radius 
of 40 pixels and subtrac ted. We fitted the X-ray spec- 
trum using the XSPEC (jArnaud 1119961 ) spectral fitting 
tool and the swxpcOtol2s6_2OO7O9OlvOll.rmf0 response 
matrix. 

Fitting the spectrum in the 0.3-10 keV band, with an 
absorbed power-law, leads to a good fit with x 2 /DoF = 
14.7/25 with a photon index of T = 1.57±g;|| and 
an equivalent hydrogen absorption column density of 
N n = 4.5± 2 ;" x 10 21 cm -2 . The absorbed flux in the 
0.3-10 keV band was 1.97xl0 -12 erg cm -2 s" 1 . The 
source clearly has a hard X-ray spectrum with a flux 
of 2.5 x 10~ 12 erg cm~ 2 s _1 , slightly lower than any 
of the previous XMM observations. The column den- 
sity i s consistent within the errors with the XMM values 
from lFeng fc KaaretJ (|2009D . The unabsorbed flux is 2.67 
xl0~ 12 erg cm -2 s _1 , corresponding to an unabsorbed 
luminosity of 4.86 x 10 39 erg s _1 in the 0.3-10 keV band 
at a distance of 3.9 Mpc. 

2.3. ATCA observations of NGC5408 X-l 

The radio nebula of NGC5408 X-l was the first de - 
tected radio counterpart of a ULX (|Kaaret et al.l 120031) . 
Later, it was con firmed that it is an extended source 
(jLang et al.l l2007j ). We obtained deep high-frequency 
ATCA CABB (Co mpact Array Broadband Backend) 
(| Wilson et al.|[201ll) observations to better constrain the 
morphology and the high-frequency part of the radio 
spectra of the nebula at 5.5, 9, and 18 GHz. 

We observed NGC5408 X-l with the CABB-upgraded 
ATCA in configuration 6D (baselines up to 6 km) on 
2009 Aug 23 (program code: C1159). The data were 
obtained simultaneously at 5.5 & 9 GHz and at 17 & 
19 GHz with 12 h on-source integration time for both sets 
(Table 1.). We observed in phase- reference mode; the 
phase calibrator was 1424-418 and the primary calibrator 
was PKS 1934-638. The data reduct ion was performe d 
using the miriad software package (jSault et al.l ll995). 
We combined the 17 GHz and 19 GHz images in order 
to enhance sensitivity. 

Fig.[5]shows the ATCA images of the nebula surround- 
ing NGC5408 X-l at 5.5, 9 and 18 GHz. The 5.5-GHz 
image was uniformly weighted, and the 9 and 18-GHz 
maps were naturally weighted in order to achieve the 
best sensitivity. Plotti ng the geometric mean beam size 
vs. the flux density, lLang et all (|2007j ) found that the 
turnover indicates that the radio emission is resolved 
with an angular size of 1.5-2.0 arcsec. Our new radio 
images show that the size of the nebula is consistent with 
the previously estimated angular size. At 9 and 18 GHz 
the source extent is only slightly larger than the corre- 
sponding beam sizes, but the elongated contours suggest 
that the nebula is resolved. This is typical for a weak, 
optically thin, steep-spectrum source, i.e. going towards 
higher frequencies, one gains resolution while the relative 
sensitivity decreases. 

We used all available measurements to fit the ra- 
dio spectrum. F ig. [S] shows previous measurements 
(jLang et alj I2007D and our new ones; combined they 
cover the 1.4-18 GHz frequency range. Our new flux 



http: / /heasarc. nasa. gov /docs/ swift / analysis / 
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Fig. 5.— The 5.5 (red), 9 (black) and 18 (blue) GHz ATCA 6D- 
array image of NGC5408 X-1. Contours represent radio emission 
at levels of 3, 4, 5, 6 and 7 times the rms noise level of 23, 17 and 
13 /jjy beam -1 , respectively. The peak brightnesses are 208, 110, 
66 /tjy beam -1 . The image was made using uniform weighting at 
5.5 GHz and natural weighting at 9 and 18 GHz and the resolutions 
are 1'.'5 x l'/O at PA=-9°, 2'.'6 x 1'.'5 at PA=-28°, l'/O x 0'.'8 at 
PA=-1° at 5.5, 9 and 18 GHz, respectively. 



NGC5408 X-1 



Frequency (Hz) 



Fig. 6.— The radio spectrum of the nebula of NGC5408 X-1. 
The best fit spectral index is a = —0.8 ±0.1 for a flux density, 

S(X!/ a . 



densities are 226±33, 137±36, and 76±20 /iJy at 5.5, 
9, and 18 GHz, respectively. The fitted radio spectral 
index of NGC5408 X-1 is a = -0.8 ± 0.1. This is the 
best constrained spectrum of a ULX nebula to date and 
is consistent with previou s results suggesting optically 
thin s ynchrotron emission (|Soria et al.ll2006aHLang et al.1 
12007ft . 

Fig. [7] shows the ATCA image at 18 GHz, overlaid on 
the HST image of t he optical nebula of NGC5408 X-1 of 
IGrise et all (|2012D . The HST filter was centered on the 
forbidden [O ill] nebular emission line (Table 2). The 
optical image shows a one-sided shell-like structure dis- 
placed from the ULX which might be due to geometrical 
effects like limb brightening towards the NE. We find that 
the radio emission at 18 GHz originates "inside" the op- 



1.0000 0.0080 0.0160 0.0240 0.0320 0.0400 00480 00560 0640 0.0720 



Fig. 7. — The 18 GHz ATCA 6D-array, naturally-weighted image 
of NGC5408 X-1 overlaid on the [O III] HST image. The sign 
'> <' marks the X-ray position of the ULX. The North direction 
is up. 

tical nebula, however it is possibly an effect of the steep 
radio spectrum. 

2.4. ESO VLT observations of NGC5408 X-1 

As one can gain information about the ionization pro- 
cess from the line flux ratios of Ha vs the forbidden 
sulphur lines, we conducted VLT observations. VLT 
FORS-2 observations of NGC5408 X-1 were obtained on 
12 April 2010 using the GRIS_1200R grism with a slit 
width of 1.0" covering the spectral range 5750-7310 A 
with dispersion 0.38 A pixel -1 and spectral resolution 
A/AA = 2140 at the central wavelengths, respectively. 
The observation block (OB) consisted of three 849 s ex- 
posures with a 12 pixel offset along the spatial axis be- 
tween successive exposures. CCD pixels were binned for 
readout by 2 in both the spatial and spectral dimensions. 
The average seeing for our new observations was 0.62 arc- 
second . Data reduction steps can be found in lCseh et al.1 
(pfTTl . Here we use a trace width of 8 pixels correspond- 
ing to 2" in order to accept a large fraction of the nebular 
emission. 

Fig. [8] shows the optical spectrum of NGC5408 X-1. 
This portion of the optical spectrum shows the forbidden 
sulphur and nitrogen lines and the Ha line. The lines are 
at wavelengths of 6558.9 ([N n]), 6573.7 (Ha), 6594.4 
QN ii]), 6727.6 ([S ii]), and 6742.0 ([S n]) in units of 
A. The corresponding fluxes are 0.46 ± 0.02, 33.1 ± 0.1, 
1.07±0.02, 1.57±0.01, 1.15±0.01 in the units of 10" 16 erg 
s -1 cm -2 . The Gaussian FWHMs of the lines represent 
the instrumental resolution of ^2.4 A. 

3. DISCUSSION 

3.1. The optical nebula around IC342 X-1 

In this section we use models developed for expand- 
ing bubbles, to calculate parameters of the nebula 
based on the observed optic al properties. We follow 
the formalism developed by iWeaver et al.1 (|1977f) and 
lOstriker &; McKeel (|1988l ) describing the hydrodynamic 
structure of a bubble. These formulae are valid for shock- 
dominated sources, i.e. jet or wind inflated bubbles, but 
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Fig. 8. — VLT optical spectrum of the nebula around NGC5408 
X-l. The dereddened flux is plotted vs. wavelength and no redshift 
correction was applied. 



not for photoionized bubbles. As IC342 X-l is considered 
to be a shock-dominated source (Pak ull fc Mirioni|[2003l: 
IRoberts et aT1l2003l; lAbolmasov et al.H2007D . we can ap- 
ply the well-known self-similar expansion law as a func- 
tion of time, t: 



R = 



125 



154tt / 



1/5 



L 



tot 

Po 



1/5 



t 



3/5 



(1) 



where L tot is the mechanical luminosity (corresponding 
to a jet or a wind or an initial explosion energy), R is the 

radius of the bubble expanding with velocity v exp = R 
into the ISM. The mass density pg is assumed to be con- 
stant and po = pm p n, where p = 1.38 is the mean atomic 
weight, m p is the proton mass and n is the hydrogen 
number density. The characteristic age of the bubble is 
t = 3R/5v exp . The kinetic energy carried by the swept- 
up mass in the expanding shell is Ek — ^L to tt, while the 
energy emitted (the cooling) by the fully radiative shock 
expanding into the ISM is E ra d — ^L to tt- The ther- 
mal energy of the gas between the reverse shock and the 



swept-up shell is E t h 



is R 



Ej^ + E r , 



li 
Et h 



Ltott- Thus, the total energy 
Ltott- 

iFeng fc Ka aret (2008) found that the bright main body 
of the nebula ha s an angular diameter o f about 6" in 
the Ha image. iPakull fc Mirionil (pOOl ): iGrise et al. I 
(12006th IFeng fe Kaaretl (|2008[) report an additional elon- 
gated structure to the South- West. Considering the en- 
tire structure of the bubble, we estimate the volume of 
the optical nebula by taking a sphere with a diameter of 
~ 10", which corresponds to ~190 pc at a distance of 
3.9 Mpc lTikhonov fc Galazutdinoval ([2010| ). 

The high flux ratio between the forbidden [S li] lines 
at 6,716 A and 6,732 A and the B ai mer Ha line indi- 
cates t he presence of shock-ionized gas. lAbolmasov et all 
(2007) infer a shock velocity of v exp ~ 20 — 100 km s _1 , 
however, line ratios of a standard library of radiative 
shock models, e.g. He n vs. H/3 is 0. 036±0.015, sugges t 
a shock velocity of ~ 100 km s" 1 (J Allen et all 12008). 
Using this velocity, we find the characteris t ic age of the 
bubble is r = 5.6 x 10 5 yr. IFeng fc Kaaretl pOOl found 



that the color-magnitude diagram suggests that the mini- 
mum stellar age in the environment of the ULX is 10 Myr. 
The characteristic age of the nebula is much shorter and 
might suggest that the nebula formation is not related 
to the formation of the central BH or BH progenitor, in- 
stead it might represent the actively accreting phase of 
the binary. 

Recalling the scaling of t he total radiative flux and 
the flux in the Balmer lines (|Dopita fc Sutherland! 1 19961 : 
lAbolmasov 6^120071^1 : 



L 



H0 



6.53 x 10" 



100 kms- 1 / 



-0.59 



(2) 



So, the total shock power represented as radiative losses 
is L ra d = 1.2 x 10 39 erg s _1 ; using the H/3 flux of 
lAbolmasov et al.l (120071 ) of 4.3 x 10- 15 erg s" 1 cm" 2 . 
As a consequence, the total mechanical luminosity is 
Ltot = ^fL rad = 3.4 x 10 39 erg s _1 , and the total kinetic 
power carried by the swept-up mass is — 6.6 x 10 38 erg 
s . Similarly, the internal (thermal) luminosity is 
L t h = 1-5 x 10 39 erg s _1 . The energy, we sec at t=r 
is E to t — LtotT — 6.0 x 10 52 erg. From Eq. 1 one can 
derive that n = 1.0 cm -3 by substituting R, r, L to t- The 
optical swept-up mass is then M = pnm p V = 2.4 x 10 38 g 
or M = 1.2 x 10 5 M s . 

IC342 X-l is some times considered as a SN remnant 
(IRoberts et all 120031 eg.). Here we intend to show that 
the total energy content does not depend significantly on 
the interpretation of the origin of the bubble. A SNR in 
the pressure-driven snowplow stage - ie. radiative domi- 
nant phase following the adiabatic phase - has an ini- 
tial explosion energy, E ACioffi. McKee fc Berts"ch inger 
H988h : 



Ei = 6.8 x 10 



43 



3.16 35 

Vkms -1 / V 



1.16 



erg 
(3) 

We have treated th e metallicity correction factor, tf; 161 , 
as 1 for clarity (|Cioffi. McKee fc Bertschingcr 1988; 
iThornton et~aH [1998) . Substituting R, v — and taking 
n = 1.0 cm~ 3 (see above), we obta in Ej = 6.0 x 10^ 2 erg. 
We note that a similar formula of iChevalierl (|1974|) pro- 
vides Ei = 5.0 x 10 52 erg. This initial energy is re- 
markably high, although somewhat expected as a sin- 
gle/simple SNR will not remain visible once it has ex- 
panded b eyond R ~ 100 p c with a canonical Ei ~ 
10 51 erg (jMatonick fc Fesenl ll997l : IRoberts et all 120031) . 
The estimated energy is in good agreement with the total 
energy content of the bubble, E to t, within model uncer- 
tainties. We note that an SNR nature for the nebula 
around IC 342 X-l is strongly challen ged by its high 
shock velocity coupled with large size ((Pakull fc Grisl 
l200l) . 

3.2. The radio nebula around IC342 X-l 

When radiation is via synchrotron emission, one can 
assume equipartition between energy of relativistic par- 



10 We note, th a t the re is a typo in Eq. 3.3 
Dopita & Sutherlan d! 111999 ) (and consequently in Eq. 1 



of 
of 



bolmasov et al. (2007)), pointed out to us by M. W. Pakull. The 
numerical factor should read 1.14 X 10 -3 rather than 2.28 X 10~ 3 
as the maximum radiative flux is ^ pv 3 an d not pv s . 
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tides and the magnetic field. We calculate the mini- 
mum total energ y of the radio nebula that corresponds 
to equipartition ([Longair II1994D : 

(4) 

where 77 — 1 is the ratio of energy in protons to rela- 
tivistic electrons, V is the volume, v is the observing 
frequency and h v is the synchrotron luminosity. As is 
customary, we do not account for relativi stic protons, 
therefore 77 - 1 = (eg. iFender et allll999Q . Substitut- 
ing the corresponding values of V ~ 1.79 x 10 56 m 3 , 
v = 5 x 10 9 Hz, L v = 3.64 x 10 18 W Hz" 1 and a fill- 
ing factor of unity, we find the energy required to power 
the nebula is E min = 9.2 x 10 50 erg. This suggests that 
the radio-emitting material carries a fraction ~ 10~ 2 of 
the initial energy. For comparison, the energy carried by 
mildly relativistic material in normal Type I c supernovae 
has b een suggested to be at most 10 -4 fPar agi et al.l 
2010), while for the jet inflated bubble around the ex- 
tragal actic microquasar S26, this fraction is a few times 
10~ 3 (ISoria et all [2010). These fractions suggest that 
most of the energy is stored in protons, nuclei, and non- 
relativistic bulk motion. 

We calculate the magnetic field s trength corresp onding 
to the minimum energy condition (jLongair 1119941) 

5 mm = 1.8xl0 10 (^) 2/ V/7 fiG (5) 

where the units and inputs are the same as above. We 
obtain B m i n — 7.4 /zG. 

When energy loss is due to synchrotron radiation, the 
lifetim e of an electron is (jLongair 1 H991 iTudose et all 
l200rl 

Substituting v = bx 10 9 Hz and B min = 7.4 //G, we find 
T sy = 18.8 Myr. So, consistently, the cooling time-scale 
is ^34 times longer than the age of the bubble and ~2 
times higher than the minimum age of the ULX stellar 
environment. 

We note that the value of rj, the energy ratio of rel- 
ativistic protons to electrons, can b e estimated from 
the high-energy pa rt of SNR spectra (jAbdo et al.ll201Ct 
lAcciari et al.l l201lT) . Assuming IC342 X-l has a ratio 
similar to that found for a scenario in which the high- 
energy gam ma-rays from the T ycho SNR are produced 
by leptons (|Acciari et al.l 120 111) , then 77 — 1 = 10 2 and 
the minimum energy obtained from Eq. [4] is increased 
by a factor of r/ 4 / 7 = 100 4 / 7 = 13.9. This would mean 
that E m i n /Etot — 0.2. Instead, if the gamma-ray emis- 
sion from Tycho is from hadrons, then r\ — 1 = 2.5 x 10 3 
and E m i n /Etot would be ~ 1.3. Thus, if the high-energy 
gamma-ray emission of shock inflated bubbles originates 
from the same electron distribution that produces the 
radio, then via the -q parameter, the value of E min could 
increase, but interestingly, does not violate the total en- 
ergy obtained from the optical. 



3.3. The radio and optical nebulae of Ho II X-l and 
NGC5408 X-l 

3.3.1. Ho II X-l 

Our new, higher resolution VLA radio observations 
clearly reveal that the morphology of the radio neb- 
ula follows the structure of the optical one (Fig. [3|). 
The radio nebula is resolved with a size of 5.5" x 
2.7", corresponding to ~ 81 pc x 40 pc. One could 
argue that this morphology r eflects either a je t activ - 
ity or an outflow. However , iPakull fc Mirionil (|2002f k 
iKaaret. Ward fc Zezasl (|2004h showed, that the nebula of 
Ho II X-l is consistent with photoionization by the cen- 
tral compact source. On the other hand, a complex ve- 
locity structure in the inner regions of the optical nebula 
indicates the im pact of the central obje ct also in the form 
of winds or jets dLehmann et al.l 12005!) . As the [S 11] vs. 
Ha ratio is <gC 1 (jAbolmasov et aLll2007f) . i.e. collisional 
excitation of the nebula is negligible, the morphology 
probably reflects an outflow rather than a well-collimated 
jet. This outflow is either relatively weaker than a jet - 
thus preventing shock-ionization - or the outflow is more 
isotropic than collimated. We note that weakly colli- 
mated outflows, in addition to jets, have been dir ectly 
observed in SS433 with VLBI fParagi et al.l 119991 ) . In 
addition, the asymmetry of the outflow is probably due 
to the fact that the nebula is den sity bounded to the East 
and South of the central obj ect (jPakull fe Mirionil 12002 1 : 
IKaaret. Ward fc Zezasl[200l . 

Furthermore, the minimum energy required to power 
the radio nebula of Ho II X-l is 2.6 x 10 49 erg with a mag- 
netic filed stren gth of 13/xG and th e synchrotron lifetime 
is r sy = 25 Myr ()Miller et al.ll2005f) . This energy require- 
ment is a factor of 35 less than needed for IC342 X-l, 
a shock-dominated source, which might also support a 
weak or uncolli mated outflow . 

Interestingly, iGrise et al. I (|2010f ) studied an X-ray 
state transition of Ho II X-l and found that it is difficult 
to interpret the thermal component of the X-ray spec- 
tra as disk emission or thermal Comptonization. On the 
other hand, this thermal component might be linked to 
the complex velocity structure of the nebula, ie. it might 
be due to a disk wind that results in a complex veloc- 
ity structure of the nearby environment. We note that 
disk wind signatures have been foun d for GRS1915 +105 
(|Neilsen fc Led 120091) and SS433 (iFabrikal l200l . i.e. 
for sources likely accreting above their Eddington-limit. 
Searching for relativistic disk lines and/or Fe K absorp- 
tion line variation in the X-ray spectrum of ULXs might 
help disentangle whether some of these ULXs are similar 
to high-accretion rate Galactic binaries. 

3.3.2. NGC5408 X-l 

Fig. E shows a shell-like optical morphology of 
NGC5408 X-l and a filled structure of the radio neb- 
ula with a maximum intensity near the ULX. The optical 
nebula is resol ved, with a diamete r of ~ 2.5" , correspond- 
ing to ~ 60 pc (jGrise et al. 1120 12T ) , which is much smaller 
than the optical nebula of IC342 X-l. Probably, the neb- 
ula is powered only by photoionization without any jet 
or outflow activity. Our optical spectra also suggest this, 
as the [S 11] vs. Ha line ratio is C 1 (Fig. [51). Fur- 
thermore, the energy needed to power the radio nebula 
of NGC5408 X-l is 3.6 x 10 49 erg with an equipartition 
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tended emission in optical and/or radio. Table |3] shows 
the main characteristics of the specific sources: the aver- 
age X-ray luminosity (Ij), the jet power estimated from 
the environment (Qj), the mass of the black hole (Mbh), 
the total energy content (E to t), the lifetime of the bubble 
(r), and the minimum energy estimated from the radio 
counterpart of the nebula (E m i n ). 

The nebula around IC10 X-l is sometimes con- 
sidered to be the result of a hypernova event 
(jLozinskava fe Moiseevl 120071 ) . However, its size is much 
larger than allowed by the surface brightness versus 
size (S-D) relation for e xplosive events, i.e. supernovae 
(jYang fe Skillmanlll993D . s uggesting that it is pow ered, 
instead, by the central BH (Bauer fc Brandtll200l . S26 
is another outlier source in the E-D relation and it 
has been revealed that the bubble is powered by a jet 
(IPakull. Soria fc MotcrJl2"o"lOl: ISoria et al.ll2010ft . There- 
fore, we will consider IC10 X-l as a black hole powered 
nebula. We note that "black hole powered nebula" means 
there is a significant and continuing outflow from the cen- 
tral black hole powering the nebula and does not argue 
against the formation of the black hole in a SN explosion. 

Considering S2 6, the X-ray photon index is T = 1.4 
(jSoria et al.ll2010f ) . consistent with a source being in the 
hard state. However, no radio core was detected with a 3- 
(7 upper limit of 0.03 mJy, which is one third of the peak 
intensity of IC342 X-l; and the total jet power of S26 is 
an order higher than for IC342 X-l (Table 3). Using the 
fundamental plane (Eq. [7]) , we obtain an upper limit on 
the mass of M ~ (1.0 ± 0.4) x 10 7 M . Given that Qj oc 
M in the hard state (iKording. Fender & MiEliarill2006bt 



field of 16jnG and a syn chrotron cooling time of r sy = 20 
Myr (|Lang et all 120071) . which is, similar to Ho II X-l, 
a factor of 25 less than the energy needed to power the 
radio nebula of IC342 X-l. 

3.4. Upper limit on the mass of the BH in ICS '42 X-l 

Self-absorbed compact jets are ubiquitous i n the 
X-ray hard states of GBH Bs ([Corbel et all l200l 
iFender. Belloni fe Gallol 12001 . IC342 X-l was found 
to have a hard X-ray spectrum in all available XMM- 
Newton and Chan dra observations covering the period 
from 2001 to 2006 (jFeng feKaaretll2009L 12008) as well as 
in our recent Swift/XRT measurement (see Sec. 2.4.1.), 
thus the presence of compact jets is expected if the hard 
X-ray spectrum equivalent to the canonical hard state of 
GBHBs. 

In addition, the morphology of the radio nebula is 
somewhat similar to the system SS433/W50 and may 
suggest a jet orientation alo ng an axis slightly East of 
North (|Feng fe Kaaretj |2b08) . As we pointed out earlier, 
an unresolved radio source is detected at the location of 
the ULX. This morphology might argue against the no- 
tion that the compact emission is the hot spot at the end 
of the jet pointing close to our line of sight. To test the 
compactness of this emission, we conducted VLBI mea- 
surements using the EVN (Sec 2.2) and we did not detect 
any source above the 3-cr noise level, thus it is likely to 
be consistent with a clump of emission from the nebula 
and might be an effect of a steep radio spectrum. We can 
set an upper limit on the flux of the putative compact 
jet of 1.9 x 10 33 erg s _1 , using our 3-cr noise level of 21 

When black holes are in the hard state, i.e. their accre- 
tion is radiatively inefficient and perhaps advection dom- 
inated or jet dominated, a relationship holds between 
X-ray luminosity, ra d io luminosity and black hole mas s 
(jMerloni et al.l &M . iFalcke.Kording fc Markoffl l200l . 
This relationship, the fundamental plane of black holes, 
has been studied on a wide mass range from black 
hole X-ray binaries to low luminosity active galactic nu- 
clei. Using the corre lation with the minimum scatter 
(jKording et al.ll2006al ): 

iog (C = L02 " ( i - 59iog (^) - iog (^) - io - 15 ) 

(7) 

and our radio and X-ray measurements, one can estimate 
the upper limit of the mass of the black hole in IC342 
X-l. However, we must consider the intrinsic rms scatter 
in the measured fundamental plane relation of 0.12 dex 
within one a. Substituting the upper limit of radio lu- 
minosity of the putative compact jet (L = 1.9 x 10 33 erg 
s _1 ) and the simultaneously measured X-ray luminosity 
Lx = 4.86 x 10 39 erg s , we estimate the mass of the 
black hole to be M BH < (1-0 ± 0.3) x 10 3 M Q . This 
limit is valid only if IC 342 X-l enters the canonical hard 
state. We note that further observations could help place 
tighter constraints on the BH mass or help test if ULXs 
exhibiting hard X-ray spectra are, indeed, in the radia- 
tively inefficient, hard X-ray state (see section 3.5.1). 

3.5. Comparison of IC342 X-l with S26, IC10 X-l 
and SS433 

In this section we compare IC342 X-l to sources whose 
jet power and accretion rate has been estimated from ex- 



iKording. Jester fc Fenderl2008l) and assuming similar ac- 
cretion efficiencies for IC342 X-l and S26; then one can 
speculate on the basis of the ratios of the jet powers 
that the average mass accretion rate of S26 is ^15 times 
higher than that of IC342 X-l. 

IC10 X-l has an X-ray photo n index of T = 1-83, po - 
tentially being in the hard state (jBauer fc Brandtl l2004) . 
Using the fundamental plane and adopting the average 
X-ray luminosity of 1.5 x 10 38 erg s _1 and a mass of ~ 30 
Mq (Table [3]), the expected core radio flux is then ^10 
/zJy at a distance of 0.7 Mpc. Thus, the source could be 
detectable with e-MERLIN or the EVLA. 

3.5.1. Jet characteristics, accretion rates and efficiencies 

In this and the following sections we investigate the 
possibility that the nebula around IC342 X-l is powered 
by a jet and the consequences regarding the jet proper- 
ties. 

The elongated morphology of the nebula and its shock- 
ionized nature are indicative of jet inflation. If the neb- 
ula is inflated by a jet, then the total power available 
in the bubble has to be equal with the time-averaged 
total jet power, ie. Qi=L to t, derived from the optica l 
bubble (IPakull. Soria fc Motchl [20101: ISoria et al.ll2010l) . 
We note that calculating the jet power using the min- 
imum energy (E m in) derived from the radio bubble 
would lead to an underestimation. This could be cither 
due to a mild deviation from equipartition, resu lting in 
E ra dio,tot/E min = 10-100 (eg. lParagi et al.ll2010l ) or due 
to the fact that the kinetic power associated to the bulk 
motio n is transferred to thermal ions also ()Soria et al.1 
[2Ol0h . 
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TABLE 3 

Comparison with X-ray sources embedded in nebula 



Name 


L x (crg/s) 


Qi (erg/s) 


M Bll (M ) 


Ltotr = Etot (erg) 


t (yr) 


E m in (erg) 


IC342 X-l 


1.6 x 10 4U 


3.4 X 10 3M 


< 1.0 x 10 3 * 


6.0 x 10 M 


5.6 x 10 b 


9.2 x 10 bu 


S26 


6.2 x 10 3b 


5 x 10 4U 


n/ a 


3.16 x 10 bi 


2 x 10 b 


10 bu 


IC10 X-l 


1.5 x 10 38 


1.27 x 10 39 


23-34 


2 x 10 52 


5 x 10 5 


n/a 


SS433 


~ 10 36 


2 x 10 38 


16 


2 x 10 51 


2 x 10 5 


10 49 ** 



Note. — The Table shows the average X-ray luminosity (Lx), the time-averaged jet power estimated from the environment (Qj), 
the mass of the black hole (Mbh), the total energy content(E tot), the lifetime of the bubble (t), and the m i nimum energy esti- 
mated from the radio counterpart of the nebulae {E m i n ). SS433 (Kirshner & Chevalier 1980; Begelman et al. 1980: Dubner ct al". 1998; 
Blundcll, Bowler & Schmidtobrcick 2008; Perez & Blundcll 2009; Fabrika 2004). S26 (Pakull, Soria & Motch 2010; Soria ct al. 2010). IC10 
X-l ( Lozinskaya & Moisccv 2007; Bauer &; Brandt 2004). *Estimated using the fundamental plane. **This value was calculated using Eq. 4 
using the parameters found bv lDubner et al.l 1 119981) : a radius of ~30 pc, a distance of 3 kpc, a radio flux density of 71 Jy at 1.4 GHz, and 
ta king a spectral index of -0.48 between 85 MHz and 5 GHz. We note that this valu e is two orders of magnitude larger than that quoted 
bv lDubner et all l|1998l ). however consistent with the value of Begelman ct al. (1980). 



In general, the total jet power is a constant frac- 
tion (/) of the ava ilable accretion power (however see 
iCoriat et al.l ()201 lh for possible variation of /). Thus, 
we can write 

Qj = fQacc = fM acc C 2 , f<l (8) 

where / is typically in the r ange of 10 ~ 3 to 10 _1 
(iFalcke fe Biermannl 11991 Il99l . 

Taking a constant rate of energy input that is charac- 
terized by the power of the jet, the rest-mass transport 
along the jet, Mj, is (jKaiser fc Alexander1ll997() : 



Taking a minimum Lorentz factor of 7 = 2 
(jMirabel fc Rodri'gued[l999l:lFendeii[2003l ). we find Mj ~ 
6. x 10~ 8 Mm yr" 1 (3.8 x 1 18 g s" 1 ). 

IKording. Jester fc Fen der (2008) found that the total 
jet power could be estimated from the flux of the compact 
jet as: 

«i - 7.2 x 10- (^^t) ergs- (10) 

This relationship was found for FR I and FR I I radio 
gala xies and scaled to C yg X-l. (See iGallo et al.l (j2005[ ) 
and iRussell et al.l (|2007t ) for the jet power estimation of 
Cyg X-l). Substituting the total jet power obtained from 
the optical nebula (Qj = 3.4 x 10 39 erg s _1 ), we find that 
the expected average luminosity of the putative radio 
core would be 6.1 x 10 33 erg s _1 . Not detecting a compact 
jet with an upper limit of 1.9 x 10 33 erg s _1 could mean 
that either its flux is variable, and currently below our 
detection limit but with an average above, or the hard 
X-ray spectrum does not represent the canonical hard 
state of GBHBs. 

In addition, it is possible to estimate the accretion 
rate, for hard state object s, based on the luminosity 
of the c ompact radio jet (Kording, Fender fc Migliaril 
(|2006bT) : IKording. Jester fc Fender! (|2008D ; and refer- 
ences therein): 

/ T \ 12 / 17 

Af ~4xl0 17 ^f^° gs -i (11) 

V10 30 ergs-V B v 1 

For IC342 X-l, we obtain an upper limit, due to the fact 
that we do not detect a compact jet, of M < 8.2 x 10 19 



g s" 1 (1.3 x 10 -6 M Q yr" 1 ). This rate is close to the 
Eddington rate for a 1OOM compact object. Comparing 
the accretion rate to the jet mass loss rate of Mj ~ 3.8 x 
1q18 g S -i leads to the conclusion that the jet efficiency 
is / > 0.046. This would be a typical value for a jet 
efficiency and could be consisten t with the assumption of 
lK5rding. Jester fc~F cndcr (2008) of f=10 _1 for obtaining 
the relationship of Eq. QTJ 

Let us compare IC342 X-l to the peculiar sources 
SS433 and GRS1915+105. The jet mass flow rate in 
SS433 is Mj = 5 x 10~ 7 M yr" 1 and the available ac- 
cretion rate is M = 10~ 4 M yr -1 (jPerez fc Blundell 
[2^TltlFabrikal[200l . so the jet efficiency is / = 5 x 10" 3 , 
which is also within the typical range for black hole jets. 
Comparing the accretion rate of IC342 X-l to SS433, we 

find Mic342X-i/-^sS433 > 10~ 2 . In contrast, the accre- 
tion rate of IC342 X-l might be an order of magnitude 
hig her than the rate of G RS 1915+105 of - 10 19 g s" 1 
(eg. iRushton et aT] 120101 ). It is interesting to note that 
if one takes the lifetime of the sur rounding bubble , and 
assumes a constant accretion rate (|Begelmanll2002T) . the 
compact object in SS433 could accrete ^20 M Q . During 
the lifetime of the nebula, the compact object in IC342 
X-l could accrete only <1 M . This comparison of 
accretion rates suggests that there is no need to invoke 
super-Eddington accretion for IC342 X-l if the mass of 
black hole is above ~100 M . 

4. CONCLUSIONS 

We studied the radio and optical nebulae of three 
ULXs. One of these ULXs, IC342 X-l, is a newly dis- 
covered radio association. 

• We confirmed that the radio spectra of the nebulae 
surrounding Ho II X-l and NGC5408 X-l are con- 
sistent with optically thin synchrotron emission. 

• We estimated the energy needed to supply both 
the optical and radio nebulae around IC342 X-l. 
The energy needed, 6 x 10 52 erg, is ~ 2 orders of 
magnitude higher than the explosion energy of a 
SN. By comparing the age of the bubble to the 
stellar environment, we found that the nebula is 
much younger, indicating that the nebula forma- 
tion is not necessarily related to the formation of 
the black hole progenitor. 



Black hole powered nebulae and the case of IC342 X-l 



11 



• We estimated that the minimum energy needed to 
power the radio nebula of IC342 X-l is 9.2 x 10 50 
erg, at least an order of magnitude higher than that 
of Ho II X-l and NGC5408 X-l. The fraction of 
energy carried by relativistic material is relatively 
high. 

• In addition to discovery of a radio nebula around 
IC342 X-l, we found a radio component unresolved 
on VLA scales, that was not detected with VLBI. 
This puts an upper limit on the flux density of a 
compact jet. The ULX was found with a hard X- 
ray spectrum in observations covering <~ 10 years. 
Use of the 'fundamental plane' of accreting black 
holes, which is valid for hard state objects, would 
place an upper limit on the mass of the black hole 
in IC342 X-l of (1.0 ± 0.3) x 10 3 M Q . 

• According to the above properties of IC342 X-l, 
we argued that the nebula is possibly inflated by a 
jet. We found that the calculated time averaged jet 
power, jet efficiency, and available accretion power 
could be consistent with that the nebula is inflated 
by a jet. 

In summary, the energetics of the surrounding nebula 
along with the possible jet properties and accretion rate 
support the idea that the nebula surrounding IC342 X-l 
could be an inflated bubble driven by the jet from the 
central black hole. 
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